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Abstract 

We have extended our earlier simulations of the high temperature behaviour 
of lattice QCD with two light flavours of staggered quarks on a 16 3 x 8 lattice 
to lower quark mass {m q = 0.00625). The transition from hadronic matter 
to a quark-gluon plasma is observed at Q/g 2 = 5.49(2) corresponding to a 
temperature of T c « 140MeV. We present measurements of observables which 
probe the nature of the quark-gluon plasma and serve to distinguish it from 
hadronic matter. Although the transition is quite abrupt, we have seen no 
indications that it is first order. 
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I. INTRODUCTION 



The understanding of the thermodynamics of hadronic matter is a major goal of 
lattice QCD. In particular, one would like to find and understand the nature of the transition 
from hadronic matter to a quark-gluon plasma and to determine the properties of this 
plasma phase which serve to distinguish it from the hadronic phase. Such high temperature 
matter must surely have existed in the early universe, and its properties would have helped 
determine the evolution of the universe. In addition, it is believed that heavy-ion collisions 
such as will be observed at RHIC will produce hot nuclear/hadronic matter and probably 
quark-gluon plasma. Finally, such studies can help us understand the dynamics of QCD 
including confinement and chiral symmetry breaking. 

We have carried out simulations of QCD with 2 flavours of light staggered quarks on 
a 16 3 x 8 lattice. In the first paper in this series we described simulations with quark mass 
m q = 0.0125, in lattice units [[|]. In the current paper we describe simulations with quark 
mass m q = 0.00625 and, where possible, attempt to extrapolate to the chiral limit. Prior 
results of this study were reported at Lattice conferences [@,[§, and there has also been a 
report from the Columbia group on their work with a lighter quark mass, 0.004 |4]]. 

The simulations were performed on the Connection Machine (CM-2) at the Pitts- 
burgh Supercomputing Center. We used the hybrid molecular dynamics algorithm with 
noisy fermions for our simulations enabling us to tune the number of quark flavours to 2. 
As at the higher mass, we found no evidence for a first-order transition, in agreement with 
work on lattices with Nt = 4 ||, 6 @, and 12 |7]]. This is what was predicted on the basis of 
chiral spin models, which predict a second-order transition at m q = and a rapid crossover, 
with no transition at small m q ||. 

The position of the transition was determined by measuring the Wilson/Polyakov line, 
which is associated with confinement, and the chiral condensate, (ipip), which measures chiral 
symmetry breaking. In addition, we have calculated the entropy densities of both quarks 
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and gluons, and measured the topological susceptibility. From the correlation functions of 
spatial and temporal Wilson lines, we have estimated the string tensions and Debye screening 
lengths, and addressed the question of ip/ J production. By measuring the baryon number 
susceptibility, we have determined how the system responds to a finite chemical potential 
for baryon number. Additionally, we have measured the hadronic screening lengths for light 
hadronic excitations. These should give us an indication as to whether the plasma has any 
hadron-like excitations or whether its excitations are simply quarks and gluons. Finally, 
we have added valence "strange" quarks, whose masses are tuned to give the correct kaon 
masses when combined with u or d quarks of mass m q = 0.00625. By measuring the entropy 
of such strange quarks, we can address the question as to whether the formation of the 
quark-gluon plasma is characterized by increased kaon production. We have also measured 
the kaon screening lengths. 

Section 2 describes the simulations and the extractions of most of the order parame- 
ters. In section 3 we discuss Debye and hadronic screening lengths. Finally, in section 4 we 
present further discussions and conclusions. 



II. SIMULATIONS 



These simulations were performed on a 16 3 x 8 lattice with 2 flavours of light quarks 
with mass m q = 0.00625, which will be compared, when appropriate, with earlier simulations 
with m q = 0.0125 [|l|. The updating was performed with time increment dt = 0.005, except 
at 6/g 2 = 6.0 and the first half of the simulations at 6/g 2 = 5.55 where it was chosen to 
be dt = 0.01. We have simulated for 1000 time units at Q/g 2 = 5.45, 1500 time units at 
6/g 2 = 5.475, 1145 time units from a cold start and 1000 time units from a hot start at 
6/g 2 = 5.5, and 1000 time units at each of 6/g 2 = 5.525, 5.55 and 6.0. Configurations were 
saved every 20 time units at 6/g 2 = 6.0, every 10 time units at 6/g 2 = 5.55, and every 5 
time units for each of the other couplings. These configurations were used for later analyses. 
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FIG. 1. The Wilson/Polyakov line (circles) and (ijjip) (crosses) as functions of 6/g 2 . 

The thermal (temporal) Wilson/Polyakov line, and the chiral condensate (ipip) were 
measured and are plotted in Fig. |I[ They are tabulated, along with values of the plaquette, 
in Table Q. In each case the first 200 time units are discarded for equilibration. The rapid 
crossover of these two quantities between 6/g 2 = 5.475 and 6/g 2 = 5.5 suggests that the 
transition occurs in this range. This is further born out by looking at the time histories of the 
Wilson/Polyakov line at 6/g 2 = 5.475 and 5.5. These are shown in Fig. At 6/g 2 = 5.5, 
starts from both the hot and cold sides of the transition evolve to the high-temperature 
phase. At 6/g 2 = 5.475, starting from the high temperature phase, the system rapidly 
settles into the low temperature phase and stays there. 

We have measured the partial entropy densities for the gluons and the light dynamical 
(u and d) quarks. In addition, we have calculated the partial entropy densities for a heavier 
"strange" (s) valence quark. At tree order the partial entropies for the quarks and gluons 
are given by 

s° 9 T=^(P st -P ss ) (1) 
3 
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TABLE I. The Wilson/Polyakov line and (ipip) as functions of 6/g 2 for m q = 0.00625 
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FIG. 2. Time evolution of the Wilson/Polyakov line (a) at 6/g 2 = 5.475 and (b) at 6/g 2 = 5.5. 
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respectively, where P st and P ss are the space-time and space-space plaquettes normalized as 
in M . Since at this order, the energy density e and pressure p obey the perfect gas relation 



p = |e, which is clearly false for a first order transition, we include the 1-loop corrections. 
Through 1-loop order ||, 
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The tree level quantities are tabulated in Table |T|. The multiple values for the strange quark 
at each f3 correspond to different choices for the strange quark mass, which span the range 
of estimates for the masses required to give the correct kaon mass when m Uj d = 0.00625. 
The values chosen were 0.035, 0.040, 0.045 and 0.050 at (3 = 5.45, 0.030, 0.035 and 0.040 at 
(3 = 5.475, 0.025, 0.030 and 0.035 at (5 = 5.5, 0.020, 0.025 and 0.030 at j3 = 5.525, 0.0175, 
0,0225, 0.0275 at (3 = 5.55 and 0.010 and 0.015 at (5 = 6.0. The 1-loop corrected partial 
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TABLE II. Partial entropy densities for gluons, u plus d quarks, and s quarks at tree order. 



entropy densities, each divided by the ratio of the corresponding free field quantity on a 
16 3 x 8 lattice to its continuum value, in an attempt to remove some of the finite size and 
finite lattice spacing effects [|TU[, are plotted in Fig. ^. The above ratios are 0.99136 for 
gluons and 1.59828 for u and d quarks. The gluon partial entropy density divided by T 4 is 
small for small values of 6/g 2 , and increases abruptly to close to the Stefan-Boltzmann value 
(32vr 2 /45) between 6/g 2 = 5.475 and 6/g 2 = 5.5, and remains there as 6/g 2 is increased. 



5 



The partial entropy for the u + d quarks over T 3 also increases rapidly between 6/g 2 = 5.475 
and 6/g 2 = 5.5, and slowly increases for higher values. Even by Q/g 2 = 6.0, it is still little 
more than half its Stefan-Boltzmann value (77r 2 n//15). This slow approach to its asymptotic 
value could well be due to poor convergence of the perturbative loop expansion. However, in 
contrast to the case of Wilson fermions, tadpole improvement is not expected to be (and in 
fact is not) of any help (Tl|. The entropy density of the strange quark behaves in a fashion 
qualitatively similar to that of the u + d quarks. Just above the transition (at 6/g 2 = 5.5) 
it is about 36% of that of the u + d quarks, while by 6/g 2 = 6.0 it has risen to about 44% 
on its way to the high temperature limit of 50%. 
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(a) (b) 
FIG. 3. Partial entropy densities of (a) the gluons and (b) the quarks as function of 6/g 2 . The 

dashed lines are the Stefan-Boltzmann values, in (b) the lower dashed line is for the s quarks and 

the upper is for the u + d quarks. 

For a complete understanding of the thermodynamics of lattice QCD one really needs 
to work at finite baryon- and hence quark-number density. Although that is beyond the scope 
of these simulations, we have measured the singlet (S) and non-singlet (NS) quark-number 
susceptibilities fl~2|| . 
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and 

XNS = {w u -Jir) {{nu) ~ {nd)) (6) 

where (n u ) and (n^) are the expectation values of the number densities for u and d quarks 
and n u and \i& are chemical potentials for u and d quarks. We have evaluated xs an d Xns 
on the lattice evaluated at zero chemical potentials using the method of \I2\ with 50 noise 



vectors (except at Q/g 2 = 5.55 where we used 20, 50, and 100 noise vectors - we report 
results only for 100 noise vectors - which indicated that 50 noise vectors was adequate), and 
plotted our results in Fig. f|. The difference between xns an d Xs remains small throughout 
the range, xs is small below the transition, rising rapidly to close to the continuum limit - 
NfT 2 = 0.03125 - as Q/g 2 is increased through the transition, exceeding the continuum value 
for large Q/g 2 , a possible finite size effect. The singlet susceptibility measures the ease with 
which a quark excess can be created by a finite chemical potential. Below the transition, 
confinement prevents the production of quarks, except as baryons, which requires extra 
energy: thus xs is small. Above the transition unconfined quarks can be created, which 
requires little energy, so xs is large, xns measures the ease with which an excess of u 
quarks and anti-d quarks can be created. Below the transition, the easiest way to do this is 
to create an excess of 7r +, s. However, the mass of the pion (larger than normal since m q is 
larger than its physical value) tends to suppress such production, and xns is small. Above 
the transition u and d's can be created with little energy and xns is large. 

Quark number susceptibility is just one of the quantities which is of importance if 
our eventual goal is to understand the equation of state near the transition. Other such 



quantities are the fluctuation quantities such as [13,14 
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FIG. 4. Quark number susceptibility xs and Xns — Xs as functions of Q/g 2 . 
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Since such quantities have peaks which span only a small range of Q/g 2 , to calculate them 
directly requires knowing their values at 6/g 2, s which are closer together than typically 
simulated. This deficiency can often be overcome by interpolating using the Ferrenberg- 
Swendsen [15 method. We found, unfortunately, that the spacing of our simulated 6/g 2, s 
- A6/g 2 = 0.025 in the transition region - was too large to permit such an interpolation. 
We estimate that one would need a A6/g 2 ~ 0.01 (or less) to enable such an interpolation 
on this size lattice (not unreasonable for future simulations). The contributions to these 
quantities at only the simulated values of Q/g 2 show little variation, and no peaks - the 
transition region is too narrow and none of our points are close enough to feel much effect. 
In the case of Xmi since only one noise vector was used for each estimate of (ipip)-, it is 
probable that the fluctuations from use of this stochastic estimator completely overwhelmed 



those associated with the gauge configurations [16 . 
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From our measurements of the Wilson/Polyakov Line, (ipip), and entropy density, 
we conclude that the transition from hadronic matter to a quark-gluon plasma occurs at 
Q/g 2 = 5.49(2) for m q = 0.00625 as compared to Q/g 2 = 5.54(2) for m q = 0.0125. Linear 
extrapolation leads us to predict that the transition would occur at Q/g 2 = 5.44(3) in the 
chiral (m q = 0) limit. This m q dependence of the position of the phase transition makes it 
difficult to extrapolate observables to m q = 0, since to extrapolate to m q = at constant 
Q/g 2 , requires Q/g 2 to lie on the same side of the transition at m q = as at m q = 0.0125. Of 
course, one might argue that it is more physical to continue, not at constant Q/g 2 but rather 
at constant physics. Since the change in Q/g 2 for the transition A(Q/g 2 ) « 0.05 between 
m q = 0.00625 and m q = 0.0125 and between m q = and m q = 0.00625, we have attempted 
linear extrapolations in mass with A(Q/g 2 ) = 0.05 between m q = 0.00625 and m q = 0.0125 
as well as for the traditional A(Q/g 2 ) = 0. In addition, we have included extrapolations for 
A (6/g 2 ) = 0.025. Such extrapolations are shown for the chiral condensate (tpip) i n Table JTf\ . 
Except at Q/g 2 = 5.45 we have restricted ourselves to those for which the m q = 0.0125, the 
m q = 0.00625 and m q = lie on the same side of the transition. At Q/g 2 = 5.45 we have 
included cases where the m q = 0.0125 and m q = 0.00625 points lie in the hadronic phase 
while the m q = point (probably) lies in the plasma phase. Although not representing the 
correct chiral limit, these, we believe, are the correct values to use to predict the topological 
charge. It is amusing to note that, above the transition (at m q = 0), where (ipip) should 
vanish at m q = 0, the extrapolations at constant Q/g 2 , even where they cross the transition, 
do better than those where Q/g 2 is changed to avoid the transition. 

We have calculated the topological charge Q on each configuration using the cooling 
method [O] described in more detail in our earlier paper M. From these Q's we have 



calculated the topological susceptibilities x as 

x = ^(Q 2 ) (io) 

The values we obtain by applying this method naively at Q/g 2 = 5.45 are actually slightly 
larger at this quark mass (m q = 0.00625) than they were at m q = 0.0125 (28(3) x 10~ 5 
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5.475 
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[0.0126(49)] 
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0.0285(24) 
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[0.0043(47)] 


0.0230(23) 
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0.0113(23) 






6.000 
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TABLE III. {ipif)} extrapolated to m q = 0. Those values marked with an asterisk are not to be 
interpreted as true m q = values, but are included because they are related to topological charge. 
Naive extrapolations which cross the transition are included in square brackets for comparison. 

compared with 24(5) x 10 -5 ) whereas x should be proportional to m q and thus smaller at 
the lower mass. For this reason we have recalculated x after removing small instantons 
by hand from the cooled configurations. Here the value at m q = 0.00625 has dropped to 
19(2) x 10~ 5 , indicating that small instantons which are a lattice artifact are part of the 
problem, which suggests that the lattice spacing is too large to allow complete separation of 



long and short distance effects on these lattices. The complete results are shown in table [TV 
and compared with the theoretical value for the confined region 
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(11) 
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which is valid for small m q . We note that, in the hadronic phase at /3 — 5.45, our measured 
value is greater than theory, as expected. In the plasma phase, we have equated theory 
to zero since (ipip) should vanish in the chiral limit in this chirally symmetric phase. The 
reason x does not, and should not vanish in this phase is that the theoretical expression 
above neglects terms higher order in m q . In fact, in the plasma phase 
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X oc m n q f = m\ (12) 

At Q/g 2 = 5.55, which is above the transition for both masses, \ = 3.0(6) for m q = 0.0125 
and, even with small instantons removed \ = 2.1(5) at m q = 0.00625, a ratio of ~ 1.4 rather 
than the predicted 4. Further analysis reveals that at least part of this is due to the finite 
lattice spacing which raises the energy of the "zero mode" associated with the instanton to 
a higher value, comparable in magnitude to the quark mass, invalidating the derivation of 
this formula |i~8| . 
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25 
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TABLE IV. Topological susceptibility x- Standard gives the results obtained by the cooling 
method. Broad are the results after removing small instantons, i.e. those where the maximum 
value of \FF\ on a site is > 2. 



We end this section with an estimate of the temperature of the transition from 
hadronic matter to a quark-gluon plasma. Once we know the critical value of Q/g 2 , the 
critical temperature T c is given by 

T < = (13) 

The lattice spacing a(g 2 ) can be determined from measurements of the hadron spectrum at 
zero temperature at 6/g 2 = Q/g 2 . Since this has not been done, we have calculated a(g 2 ), 
from a value of the p mass obtained by interpolation from simulations at other values of 
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6/g 2 and m q , following Blum et al. 0. This yields T c = 150(9)MeV for m q = 0.0125, 
T c = 140(8)MeV for m g = 0.00625, and T c = 128(9) at m q = 0. 



III. SCREENING LENGTHS 

In this section we present our measurements of screening masses, the inverse of screen- 
ing lengths. We measure the correlations of time-like Wilson lines which are related to the 
string tension in hadronic matter, and the Debye screening length/mass in the quark-gluon 
plasma, which is the range of influence of a static colour charge in the plasma. The corre- 
lations of spacelike Wilson lines are also measured. In addition we measure the screening 
masses for excitations with the quantum numbers of low mass hadrons, to determine whether 
quark-gluon plasma has hadron-like excitations, or whether it is best described in terms of 
quark and gluon excitations. 

As described in , we measure correlations between "fuzzy" Wilson/Polyakov loops 
with zero x, y, and t momenta, correlated in the z direction. Again, we find that those which 
were "blocked" the maximum number of times (3), i.e. maximally "fuzzy" were the most 
correlated. For timelike Wilson lines our subtracted correlation function should behave at 
large z separations Z as 

P(Z) = A{exp[-fiZ} + exp[-/2(N s - Z)}}. (14) 

As usual one obtains effective screening masses n(Z) (where we have suppressed the mass 
dependence) by fitting this form for PiZ — 1) and P{Z). Then fi(Z) — * /i, the screening 
mass, which is closely related to the Debye screening mass PD|| . The screening masses so 
obtained are tabulated in Table [V]. As can be seen, except for the 6/g 2 = 5.55 results, our 
data show some evidence for a plateau as early as the Z=l-2 effective masses (/i(2)), so we 
use /i(2) as our estimate for /i. In the plasma phase, the quantity of interest /x/T has the 
values 3.8(7) at 6/g 2 = 5.5, 4.8(5) at 6/g 2 = 5.525, 6.0(1.2) at 6/g 2 = 5.55 and 5.4(5) at 
6/g 2 = 6.0, not significantly different from the m q = 0.0125 results. 
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Temporal Wilson/Polyakov Line Effective Mass 


6/g 2 


Z=0-1 


Z=l-2 


Z=2~3 


Z=3-4 


5.450 


1.10(8) 


0.84(11) 


0.96(26) 


1.35(80) 


5.475 


0.94(4) 


0.70(6) 


0.73(15) 


0.93(46) 


5.500 


0.70(8) 


0.49(9) 


0.49(13) 


0.57(30) 


5.525 


0.77(5) 


0.60(6) 


0.66(13) 


0.86(41) 


5.550 


0.84(7) 


0.75(15) 


0.91(40) 




6.000 


0.74(4) 


0.68(6) 


0.66(16) 


0.78(41) 



TABLE V. Effective masses fJ,(Z) from temporal Wilson line correlations. 



The correlations C(Z) between spatially oriented "fuzzy" Wilson/Polyakov lines, also 
have an exponential asymptotic behaviour. Unlike the temporal Wilson/Polyakov lines, the 
spatial lines do not develop a spatial expectation value corresponding to a perimeter law 
in the high temperature phase, but retain an area law defining a confining theory in both 
phases. The effective masses we obtained from C(Z) for maximal blocking are given in 
table [VI]. We note that these values are very similar to those at m g = 0.0125. Taking /i(2), 







Spatial Wilson/Polyakov Line Effective Mass 


6/g 2 


Z=0-1 


Z=l-2 


Z=2-3 


Z=3~4 


5.450 


2.80(7) 


1.93(43) 






5.475 


2.45(6) 


1.68(22) 


0.91(83) 




5.500 


1.92(5) 


1.10(12) 






5.525 


1.74(4) 


1.19(12) 


1.28(34) 




5.550 


1.53(5) 


0.97(8) 


1.09(31) 




6.000 


0.51(2) 


0.41(3) 


0.39(3) 


0.40(5) 



TABLE VI. Effective masses from spatial Wilson/Polyakov line correlation functions 



the value from fitting C(Z) from Z = 1 to Z = 2, as our best estimate of the mass //, we 
can extract a string tension k = fx/N t . Fig. |] shows y^/m^ as a function of 6/g 2 , indicating 
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that the string tension appears to be constant in physical units. This is in accord with 
expectations that these spatial Wilson lines should show confinement at all temperatures 
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FIG. 5. Spatial string tension function of 6/g . 

Let us now turn to the hadron screening lengths, or rather their inverses, the hadron 
screening masses. Not only can measurements of these screening masses give indications of 
the nature of excitations in hadronic matter and the quark-gluon plasma, but they can also 
shed light on the nature of chiral symmetry restoration JH]]. A full understanding of chiral 
symmetry realizations at finite temperature requires a study of the flavour singlet mesons 
with their disconnected contributions, which is beyond the scope of this paper. However, as 
we will see, we can at least address some of these issues. 

For massless quarks, the SU(2)a component of the SU(2) x SU(2) chiral symmetry 
is spontaneously broken, while the SU(2)v symmetry is unbroken. The U(l)v symmetry 
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(baryon/quark number) is conserved, while the U(1)a is broken explicitly by the chiral 
anomaly, as well as spontaneously At the finite temperature transition, the chiral condensate 
vanishes, so one expects the SU(2)a and hence the 577(2) x SU(2) chiral symmetry to be 
restored. What is less clear, is whether the U(1)a is restored to give U(2) x U(2) chiral 
symmetry 



The description of which spatial hadron propagators were considered to get our 
screening masses was given in our first paper, and will not be repeated here, except to 
note that wall sources and the x — y — t Coulomb gauge were employed again, the correla- 
tions being measured in the z direction. We also measured screening masses for hadronic 
excitations involving valence "strange" quarks whose masses are given above. 

Fig. |6]a shows the screening masses of the 7r and the connected part of the fo (cr) 
propagator - the do propagator. The error bars reflect only statistical errors. These two 
masses would become equal when U(1)a symmetry is restored. The ir and ao appear to 
have an appreciable mass splitting just above the transition, only approaching one another 
for considerably higher 6/g 2 . The fact that this splitting appears to be greater than at 
m q = 0.0125 indicates that finite quark mass is probably not the reason for this splitting. 
However, since the ao mass is difficult to measure we cannot rule out at least part of this 
splitting being due to systematic errors. More importantly, since flavour symmetry breaking 
is large at these 6/g 2 values, we cannot rule out the splitting being due to flavour symmetry 
breaking. Taken at face value this data suggests that U(1)a symmetry is not restored at 
the transition, but instead is restored as T — > oo. Recent measurements of the disconnected 
part of the fo propagator QTSj , and the susceptibility associated with it p3|-^5|] indicate 



that it splits the fo and a masses in the neighbourhood of the transition and gives fo- 
il degeneracy. This supports the scenario where SU{2)a is restored at the transition, as 
expected, but U(1)a is only restored at somewhat higher temperatures. 

The p and a± screening masses as a function of 6/g 2 are presented in Fig. |6|b. Since 
these states lie in the same SU(2) x SU{2) multiplet, their masses should become equal 
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FIG. 6. Hadronic screening masses as functions of 6/g 2 : (a) Scalar states - tt, ao, (b) Vector 
states - p, oi, (c) Nucleon states - N, N', (d) ET, tt. 

when iS f Z7"(2) > 4 is restored. What we observe is that these screening masses do come together 
just above the transition indicating that SU(2) x SU(2) is restored at the transition, in the 
chiral limit. Similarly, Fig. ^c indicates that the N and N' (negative parity nucleon) become 
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degenerate at (or near) the transition which is also evidence that chiral SU(2) x SU(2) 
symmetry is restored. 

We also note that, just above the transition, the p and a\ masses lie just below the 
energy of two free quarks at the lowest Matsubara frequency, i.e. 2ttT, and presumably 
approach this value from below as T — > oo. Similarly the N and N' masses lie just below 
the energy of three free quarks at the lowest Matsubara frequency (3irT), and approach 
this value at high temperatures. Thus these states appear to be weakly bound above T c . In 
contrast, both the ir and the masses are considerably below 2nT, just above the transition, 
only approaching this limit as T — > oo, indicating that they remain strongly bound, even 
in what is usually thought of as the deconfmed phase. Thus the quark-gluon plasma does 
appear to support hadronic excitations. 

Finally in Fig. ^|d, we present our kaon screening lengths for the choices of the valence 
strange quark mass given above. The pion screening lengths are given for comparison. Note 
that, the pion screening mass is for m u ^ = 0.00625 which is considerably larger than the 
physical quark mass, so that the pion screening masses will be too large. However, the 
strange quark masses are tuned to give reasonable values for the kaon mass, so that kaon 
screening lengths we obtained should be reasonable estimates for the physical kaon screening 
lengths. 

IV. DISCUSSION AND CONCLUSIONS 

We have extended our previous simulations of the thermodynamics of lattice QCD 
with two light staggered quarks down to quark mass m q = 0.00625 (lattice units), in an 
attempt to probe the chiral limit. At this mass we find a transition at a critical coupling 
6/<7g = 5.49(2) corresponding to a critical temperature T c = 140(8)MeV. Combining this 
with our previous result at m q = 0.0125, we predict a critical coupling given by Q/gl = 
5.44(3) corresponding to T c = 128(9)MeV in the chiral (m q = 0) limit. As expected, there 
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is no evidence for a first order transition. 

The screening length associated with thermal Wilson/Polyakov line correlations is 
quantitatively similar to that obtained from our published m q = 0.0125 calculations. Thus 
our conclusion that we should expect ip/J suppression for T > 2T C still stands, while the 
situation for T c < T < 2T C remains unclear. 

We observed that the entropy of strange quarks is ~ 70% of that of either u or d 
quarks, just above the transition. It is to be hoped that phenomenology based of this number 
might shed light on whether there is an increase in kaon production as the system passes 
through this transition. 

Our study of hadron screening lengths indicates that the ir, K and ao states remain 
strongly bound above the transition. Thus there are excitations with these quantum numbers 
(and presumably those of the / ), in the quark-gluon plasma. Such binding is presumably 
related to the finite string tension we have observed in our spatial Wilson line correlations. 
The other hadronic excitations we observed (p, a±, N and N') appear to be only weakly 
bound above the transition. The relatively large separation of the n and ao masses above the 
transition compared with that of the p and a x suggests that these do not become degenerate 
at the transition, so that the restored symmetry above the transition is SU (2) x SU (2) x 
LV(1) rather than U(2) x U(2) as some have suggested. This is in agreement with recent 
calculations of the chiral susceptibility Xm, and calculations of the fo mass from the full 
(connected + disconnected) propagator. 

Our calculation of the baryon number susceptibility, gives us a glimpse of physics 
at finite baryon number density. Calculation of topological susceptibility at this low mass 
appears to be strongly affected by the large lattice spacing at T c , preventing us from reliably 
extracting its mass dependence. 

For the standard Wilson action which we used, we need a larger value of N t to 
remove finite lattice spacing effects. We also need a larger aspect ratio (N s /N t ) to make 
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finite lattice size correction factors closer to one, reducing the systematic uncertainties in 
extracting entropy densities, energy densities and pressures. The larger N s is also needed 
to reduce the systematic errors involved in extracting screening lengths. Zero temperature 
values of the plaquette observable and (ipip) are needed to extract partial pressures and 
energy densities. Measurements closer together in 6/g 2 will be necessary if we are to extract 
the fluctuation quantities involved in obtaining the equation of state. 
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